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Abstract 

I review the subjects of non-solar cosmic rays (CRs) and long-duration gamma- 
ray bursts (GRBs). Of the various interpretations of these phenomena, the one 
best supported by the data is the following. Accreting compact objects, such as 
black holes, are seen to emit relativistic puffs of plasma: 'cannonballs' (CBs). 
The inner domain of a rotating star whose core has collapsed resembles such an 
accreting system. This suggests that core-collapse supernovae (SNe) emit CBs, 
as SN1987A did. The fate of a CB as it exits a SN and travels in space can be 
studied as a function of the CB's mass and energy, and of 'ambient' properties: 
the encountered matter- and light- distributions, the composition of the former, 
and the location of intelligent observers. The latter may conclude that the 
interactions of CBs with ambient matter and light generate CRs and GRBs, all 
of whose properties can be described by this 'CB model' with few parameters 
and simple physics. GRB data are still being taken in unscrutinizcd domains 
of energy and timing. They agree accurately with the model's predictions. CR 
data are centenary. Their precision will improve, but new striking predictions 
are unlikely. Yet, a one- free-parameter description of all CR data works very 
well. This is a bit as if one discovered QED today and only needed to fit a. 



1 Introduction 

This is a version of an introductory talk to high-energy physicists. Cosmic 
rays (CRs) were the first item in their field, and will remain the energy record- 
breakers for the foreseeable future. I shall argue that nothing 'besides the 
standard model' is required to understand CRs of any energy, subtracting from 
their interest. 'Long' gamma-ray bursts (GRBs) are flashes of mainly sub- 
MeV photons, originating in supernova (SN) explosions. The 7-rays are highly 
collimated. Hence, GRBs are not the publicized 'highest-energy explosions after 
the big bang', but more modest torches occasionally pointing to the observer. 
GRBs arc of interest because their understanding is intimately related to that 
of CRs. It might have been more precise to say 'my understanding' of GRBs 
and CRs, for the work of my coauthors and I is viewed as unorthodox. 

What a start! I have already admitted that our stance is not trendy and 
that the subject is of no post-standard interest. But our claims arc based on 
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Figure 1: Left (a) The weighed all-particle CR spectrum E 3 dF/dE. Some of 
the high-energy data disagree with others. Right (b) The low-energy H an He 
flux versus the kinetic energy per nucleon, multiplied by E^. In (a) the lines 
are two extreme CB-model predictions. In (b) both lines are predictions. The 
data coincide best with the prediction at solar minimum. The normalizations, 
shown here to coincide with the data, are predicted to within a factor of ~3. 



clear hypothesis, which may be proven wrong, and very basic physics, which is 
precise enough, very pretty, understandable to undergraduates, and successful. 

The information about GRBs and CRs is overwhelming. GRBs are known 
since the late 60's and CRs since 1912. Surprisingly, no theories have arisen 
that are both accepted ('standard') and acceptable (transparent, predictive 
and successful). I cannot refer to a representative subset of the ^70 + 70 kilo- 
papers on CRs and GRBs. For reviews of the standard views on CRs, see 

e.g. Hillas © or Hoerandel © . For the accepted 'fireball' model of GRBs see 
e.g. Meszaros © or Piran©. Fewer self-citations and many more references, 
particularly to data, appear in DDD02©, DDD03®, DDCM^and DD06©. 

2 Most of what you may want to recall about Cosmic Rays 

In CR physics 'all-particle' refers to nuclei: all charged CRs but electrons. 
The CR spectra being fairly featureless, it is customary to weigh them with 
powers of energy, to over-emphasize their features. The E 3 dF/dE all-particle 
spectrum is shown in fig. [T£i, not updated for recent data at the high-energy 
tail. At less than TeV energies the CR flux is larger than 1 m~ 2 s _1 sr _1 and it 
is possible to measure the charge Z and mass number A of individual particles 
with, e.g., a magnetic spectrometer in a balloon, or in orbit. Some low-energy 
results for H and He are shown in fig. [TJj. They vary with solar activity. 
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Figure 2: The CR fluxes at 1 TeV and the relative abundances in the ISM. 
Left (a) Fluxes of elements up to Z and abundances in the solar ISM. Right 
(b) Solar ISM abundances (stars), CR fluxes (squares) of the primary CRs, 
and the corresponding CB- model predictions (circles), all normalized to H. 



The CR fluxes of the lightest 30 elements at E= 1 TeV (of a nucleus, not 
per-nucleon) are shown in fig. [5^, and compared with the relative abundances 
in the interstellar medium (ISM) of the solar neighbourhood. Elements such 
as Li, Be and B are relatively enhanced in the CRs, they result from colli- 
sional fragmentation of heavier and abundant 'primaries' such as C, N and O. 
Otherwise, the solar-ISM and CR Z-distributions are akin, but for H and He. 
In fig. (2}d the abundances relative to H of the CR primary elements up to Ni 
(Z = 28) are plotted as blue squares. The stars are solar ISM abundances. CR 
positrons and antiprotons attract attention as putative dark matter products, 
but it is nearly impossible to prove that their fluxes are not entirely secondary. 

The Galaxy has a complex magnetic-field structure with B G = 0(1 fiG) 
and coherent domains ranging in size up to ~ 1 kpc, ~ 1/8 of our distance to 
the Galactic center. In such a field, a nucleus of E~p> Z (3x 10 9 ) GeV would 
hardly be deflected. For Z = 1, this energy happens to be the 'ankle' energy, 
at which the flux of fig. [1] bends up. CRs originating within the Galaxy and 
having E>E all k\ c would escape practically unhindered. The CR flux does not 
bend down at that energy, thus the generally agreed conclusion that CRs above 
the ankle are mainly extragalactic. CRs of Galactic origin and E < -E a nkie are 
'confined', implying that their observed and source fluxes obey: 

fo«r con ,4 r conf cx {Z/pY°°» f , /3 conf ~0.6±0.1, (1) 

with T C onf a 'confinement time', deduced from the study of stable and unstable 
CRs and their fragments. 

At E = 10 6 -10 8 GeV the all-particle spectrum of fig. [T^ bends in one or 



two 'knees'. The knee flux is too small to measure directly its energy and com- 
position, which arc inferred from the properties of the CR shower of hadrons, 
7's, e's and fi's, initiated by the CR in the upper atmosphere. The results for 
H, He and Fe are shown in fig. [3J Note that even the same data leads to in- 
compatible results, depending on the Monte-Carlo program used to analize the 
showers. But the spectra of the various elements seem to have 'knees' which 
scale roughly with A or Z, the data not been good enough to distinguish. 

The high-i? end of the i? 3 -weighed CR spectrum is shown in fig. 
These data and the more recent ones of HIRES and Auger, clearly show a 
cutoff, predicted by Greisen, Zatsepin and Kuzmin (GZK) as the result of 
the inevitable interactions of extragalactic CR protons with the microwave 
background radiation. The reactions p + "f — > n + n + ; p + ir° cut off the flux at 
E > E GZK ~ix 10 11 GcV, from distances larger than tens of Mpc. Similarly, 
extragalactic nuclei of E > 10 9 GeV are efficiently photo-dissociated in the 
cosmic infrared radiation, the corresponding CR flux should not contain many. 

At very high energies, rough measures of the CR A-distribution are ex- 
tracted from the 'depth of shower maximum', X, the number of grams/cm 2 of 
atmosphere travelled by a CR shower before its e ± /7 constituency reaches a 
maximum. At a fixed energy, X decreases with A, since a nucleus is an easily 
broken bag of nucleons of energy ~ E/A. As in fig. |4)d, the data are often 
presented as (ln[A](E)) , which approximately satisfies X(A)~ X(l) — xln[A], 
with cc~37 grams/cm 2 the radiation length in air. 

If CRs are chiefly Galactic in origin, their accelerators must compensate 
for the escape of CRs from the Galaxy, to sustain the observed CR flux: it is 
known from meteorite records that the flux has been steady for the past few 
Giga-years. The Milky Way's luminosity in CRs must therefore satisfy: 

L = — f — E^- dEdV ~ 1.5 xlO 41 ergs" 1 , (2) 
c J W dE 

where V is a CR-confinement volume. The quoted standard estimate of L cu is 
very model-dependent ^3 . 

3 More than you ever wanted to know about Gamma-Ray Bursts 

Two 7-ray count rates of GRBs, peaking at dN/dt = O(l0 4 ) s _1 , are shown in 
fig-GD The typical energy of the 7-ray of GBBs is ~ 250 keV. The total 'isotropic 
equivalent' energy of a source of such photons at a typical redshift, z = 0(l), is 
£^ SO ~10 53 erg, similar to the available energy in a core-collapse SN explosion, 
i.e. half of the binding energy of a solar-mass neutron star, maybe a bit more 
for a black-hole remnant. It is hard to imagine a process with > 1% efficiency 
for 7-ray production. Since GRBs are observed to be made by SNe, either the 
parent stars are amazingly special, or the 7-rays are narrowly beamed. 
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Figure 3: E 2 ^ dF/dE CR spectra in the 'knee' region. Top: protons; middle: 
a particles; bottom: iron nuclei. The data ensemble was kindly provided by 
K.H. Kampcrt. The CB-model predictions are explained in the text. Notice 
the Fe 'hip', occurring at the same 7 as the H, He hips of fig. [TJj. 



Figure 4: Left (a) The i? 3 -weighed CR flux at the highest observed energies. 
The data are not the oldest, nor the most recent, but show the currently ob- 
served trend. Right (b) The CR mean of the Neperian log of A, versus energy. 
In (a) the lines are two extreme CB-model predictions. The normalizations be- 
low and above the ankle, shown here to coincide with the data, are separately 
predicted to within a factor of ~3. In (b) the lines are as in (a). 

The total-duration distribution of the 7-rays of GRBs has two peaks, 
with a trough at ~2s dividing (by definition) two distinct types. 'Long' GRBs 
are more common and better measured than short ones; one is more confident 
discussing mainly the former, as I shall. The long GRB light curves of fig. [5] are 
not atypical. The 'pulses' of a given GRB vary in intensity, but have similar 
widths, a fairly universal exponential rise, and a power decay oct~ a , a ~2. The 
number of 'clear pulses' averages to ~5, it may reach ~12. The pulse-to-pulse 
delays are random, extending from 0(ls) to O(10 2 s). Put all the above in a 
random-generator and, concerning long GRBs, 'you have seen them all'. 

GRBs are not often seen more than once a day, they are baptized with 
their observation date. GRBs 980425 and 030329, shown in fig.0 originated at 
z = 0.0085 (the record smallest) and 0.168, respectively. How are the redshifts 
known? GRBs have "afterglows" (AGs) : they are observable in radio to X-rays 
for months after their 7-ray signal peters out. The AG of GRB 030329 in the 
'R-band' (a red-light interval) and radio is shown in fig. [S^-c. Once the object 
is seen in optical or radio, its direction can be determined with much greater 
precision than via 7 rays. Very often the source is localized within a galaxy, 
whose lines can be measured to determine z (in some cases a lower limit on z 
is deduced from absorption lines in intervening material). 

GRB 980425 was 'associated' with a supernova called SN1998bw: within 
directional errors and within a timing uncertainty of <~ 1 day, they coincided. 
The luminosity of a 1998bw-likc SN peaks at ~ 15 (1 + z) days. The SN light 
competes at that time and frequency with the AG of its GRB, and it is not 



"a Trigger 2Q030329_1 1 371 4.70 = 131558670.9 = 732973047.70 
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Figure 5: The 7-ray count-rate light curves of GRB 980425 (a: left) and GRB 
030329 (b: right). In the CB model, each pulse corresponds to one cannonball. 
The single pulse in (a) and the two pulses in (b) are fit with Eq. |H 



always easily detectable. Iff one has a predictive theory of AGs, one may test 
whether GRBs are associated with 'standard torch' SNe, akin to SN1998bw, 
'transported' to the GRBs' redshifts. The test was already conclusive (to us) 

111 2001 © . One could even foretell the date in which a GRB's SN would be 
discovered. For example, GRB 030329 was so 'very near' at z = 0.168, that one 
could not resist posting such a daring prediction during the first few days 
of AG observations. The prediction and the subsequent SN signal are shown in 
fig-Ekib. The spectrum of this SN was very well measured and seen to coincide 
snugly with that of SN1998bw, and this is why the SN/GRB association ceased 
to be doubted: long GRBs are made by core-collapse SNe. 

Astrophysicists classify SNe in Types, mainly depending on the compo- 
sition of their ejecta. Within very limited statistics the SNe associated with 
GRBs are of Type Ib/c. These constitute some 15% of core-collapse SNe, the 
fascinating ones which beget neutrinos, neutron stars and presumably black 
holes. Type la SNe are probably mere explosions of accreting white dwarfs, 
but they are very luminous, and of cosmological standard-candle fame. 

GRBs have many 'typical' properties. Their spherical-equivalent number 
of 7-rays is ~10 59 . Their spectrum at fixed t is very well approximated by: 
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with 6~1, a~l, /3~2.1. Early in the evolution of a pulse, the 'peak energy' 
(characterizing the photons carrying most of the GRB's energy) is E p ~T[0]~ 
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Figure 6: Left to right and top to bottom, a) and b) The R-band AG of GRB 
030329. A micro- Jansky is 10~ 29 erg cm -2 s -1 . a) Six days of data are used 
to predict the next ~ 100 days, and the SN contribution, b) The SN is seen, c) 
The radio afterglows of GRB 030329. d) The predicted and observed inter-CB 
superluminal angle in this GRB. (e,f): SN1987A and its two 'mystery spot' 
CBsE}. xhe motion of the Northern one was superluminal. 



250 keV, evolving later to T(t)~i 2 . A pulse's shape at fixed E is well fit by: 
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with A.t~0(\s) at E^E V . Eq. ^ reflects an approximate spectro-temporal 
correlation whereby E dN/(dE dt)~ F[Et 2 ], which we call the Et 2 'law'. 

The values of E p ; of the isotropic-equivalent energy and luminosity, E^° 
and Lp°; of a pulse's rise-time t r - lse ; or of its 'lag-time' ti ag (a measure of how 
a pulse peaks at a later time in a lower energy interval) vary from GRB to 
GRB over orders of magnitude. But they are strongly correlated, as shown 
in Figs. ([7k.-d). It is patently obvious that such an organized set of results is 
carrying a strong and simple message, which we shall decipher. 

X-ray flashes (XRFs) are lower-energy kinsfolk of GRBs. They are defined 
by having E p < 50 keV. Their pulses are wider than the ones of GRBs and their 
overlap is more pronounced, since the total durations of (multi-pulse) XRFs 
and GRBs are not significantly different. In fig. [7J; I show the time at which the 
single pulse of XRF 060218 peaked (measured from the start of the count-rate 
rise) as measured in different energy intervals. This is an impressive validation 
of the Et 2 law (the red line), also screaming for a simple explanation. 

Analytical expressions summarizing the behaviour of GRB and XRF af- 
terglows in time (from seconds to months) and frequency (from radio to X-rays) 
do exist (DDD02/03), but they are somewhat more complex than Eqs. (|3I4[) . 
The typical AG behaviour is shown in fig. [7f, as a function of frequency, at 
1, 10... 300 days after burst (the value of p is ~ 2.2 ± 0.2). This simple fig- 
ure reflects a rich behaviour in time and frequency. 'Chromatic bends' (called 
'breaks' in the literature) are an example. At a fixed time, the spectra steepen 
from ~ ^~ - 5 to ~ v~ 1A at the dots in the figure. Around a given frequency, 
such as the optical one marked by a dotted line, the optical spectrum makes 
this same transition as a function of time (at i~3d, for the parameters of this 
example), while the spectral shape at X-ray frequencies stays put. 

3.1 The Swift era 

Physicists, unlike ordinary year-counting mortals, live in 'eras'. Many are wait- 
ing for the LHC era or the Plank era, GRB astronomers are in their 'Swift era'. 
Various satellites are currently contributing to a wealth of new data on GRBs 
and XRFs. Swift is one of them. Within 15 seconds after detection, its 15-150 
keV Burst Alert Telescope sends to ground a 1 to 4 arcmin position estimate, 
for use by robotic optical ground telescopes. In 20 to 75 s, Swift slews to bring 
the burst location into the field of view of its 0.3-10 keV X-ray Telescope and 
its 170-650 nm UV/Optical Telescope. With nominal celerity, Swift has filled 
a gap in GRB data: the very 'prompt' X-ray and optical radiations. 





Figure 7: Left to right and top to bottom: a) The [E^°,Lp°] correlation, b) 
The [E p , El^°] correlation. The limiting slopes are in red, the central predictions 
in blue, c) The [t lag7 L^°] correlation, d) The [t™,Lj, so ] correlation. The data 
are for Swift-era GRBs of known z. e) The peak times of the one-pulse X-ray 
flash XRF 060218 at different energy intervals, f) A typical AG spectrum at 
different number of days after the burst (1, 10,... 300). In the figure p/2~ 1.1. 



Swift has established a canonical behaviour of the X-ray and optical AGs 
of a large fraction of GRBs. The X-ray fluence decreases very fast from a 
'prompt' maximum. It subsequently turns into a 'plateau'. After a time of 
O(ld), the fluence bends (has a 'break', in the usual parlance) and steepens 
to a power-decline. In fig. [5£i, this is shown for a Swift GRB. This bend is 
achromatic: the UV and optical light curves vary in proportion to it. Although 
all this is considered a surprise, it is not. In fig.EJa I show a pre-Swift AG and 
its interpretation in two models. In fig. [5J: one can see that the bend of this 
GRB was achromatic. Even the good old GRB 980425, the first to be clearly 
associated with a SN, sketched a canonical X-ray light curve, see fig. [5J1. 

The 7 rays of a GRB occur in a series of pulses, 1 and 2 in the examples of 
fig. Swift has clearly established that somewhat wider X-ray flares coincide 
with the 7 pulses, having, within errors, the same start-up time. On occasion, 
even wider optical humps are seen, as in fig. The X-ray counterpart of the 
second hump in this figure is clearly seen in fig. [5Jd. In an XRF the X-ray 
flares can be very wide, as in the one-flare example of fig. [TUb . In such a case, 
the accompanying optical 'humps' peak very late, at t = 0(ld), as in fig. [TUb. 
All these interconnected 7-pulses, X-ray flares and optical humps are described 
by Eqs. (|3l4p . They are obviously manifestations of a common underlying 
phenomenon, which we shall dig out. Finally, Swift has discovered that not all 
X-ray light curves are smooth after the onset of their fast decay, as the one in 
fig- HUk is. Well after 7 pulses are no longer seen, relatively weak X-ray flares 
may still be observable, as is the case in figs. ISbfTUH. 

4 Breath-taking entities: the astrophysical jets of cannonballs 

A look at the web -or at the sky, if you have the means- results in the realization 
that jets are emitted by many astrophysical systems. One impressive case 
is the quasar Pictor A, shown in figs. (fTTb .bV Somehow, its active galactic 
nucleus is discontinuously spitting something that does not appear to expand 
sideways before it stops and blows up, having by then travelled almost 10 6 light 
years. Many such systems have been observed. They are very relativistic: 
the Lorentz factors (LFs) j = E/(mc 2 ) of their ejecta are typically of O(10). 
The mechanism responsible for these mighty ejections — suspected to be due to 
episodes of violent accretion into a very massive black hole — is not understood. 

In our galaxy there are 'micro-quasars', whose central black hole's mass is 
a few Mq. The best studied is GRS 1915+105. In a non-periodic manner, about 
once a month, this object emits two oppositely directed cannonballs, travelling 
at v ~ 0.92 c. When this happens, the continuous X-ray emissions — attributed 
to an unstable accretion disk — temporarily decrease. Another example is the 
/i-quasar XTE J1550-564, shown in fig. [TTb . The process reminds one of the 
blobs thrown up as the water closes into the 'hole' made by a stone dropped onto 
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Figure 8: Right to left and top to bottom, a) The canonical 0.2-10 keV X-ray 
light curve of GRB 050315, fit to the CB model, b) Pre-SWIFT predictions for 

the 2-10 keV X-ray AG in the CB (DDD02) and fireball^ models, compared 
to data for GRB 990123. c) Broad band optical data on GRB 990123, fit in the 
CB model (DDD03). The 'bend' in (b) and (c) is due to the CBs' deceleration, 
and is achromatic, d) The X-ray light curve of GRB 980425 showing a very 
pronounced 'canonical' behaviour and what we called (DDD02) a long 'plateau'. 
The last (red) point postdates the original figure. 




Time since trigger (s) Time since trigger (s) 

Figure 9: Left to right and top to bottom: a) Early R-band light curves of GRB 
060206. b) Its X-ray light curve, c) The X-ray light curve of GRB 060904a 
with its late flares: the progressively dying pangs of its accreting engine, d) 
The hardness ratio of GRB 060904a traces the 'prompt' ICS pulses of its light 
curve, settling to a constant as SR becomes dominant in the 'afterglow'. The 
above understanding of all these data is specific to the CB model. 
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Figure 10: Left to right and top to bottom: a) The X-ray light curve of XRF 
060218. b) Data on XRF 060218/SN2006aj. Upper part: the 0.3-10 keV 

SWIFT-XRT light curve, with fits by Campana et al. Q3>. Lower part: UVO 
light curves. In our model, the X-ray flare and optical humps are made by ICS 
by a single CB. c) The hardness ratio of XRF 060218. d) The extensive X-ray 
light curve of GRB 061121. The lines in (a,c,d) are CB-model fits. 




Figure 11: Left to right and top to bottom, a) and b) X-ray images of Pictor 
A. b) Image centred at the leftmost spot in (a) and superimposed on VLA 
radio contours, c) Two relativistic CBs emitted in opposite directions by the 
microquasar XTE J1550-564, seen in X-rays, d) HST picture from 28 October 
2002 of the glory, or light echo, of a pre-supernova outburst of the red supergiant 
V838 Monocerotis, doctored with some CB-model art-work. 



its surface, but it is not understood; for quasars and /x-quasars, the 'cannon's' 
relativistic, general-relativistic, catastrophic, magneto-hydro-dynamic details 
remain to be filled in! Atomic lines of many elements have been seen in the 
CBs of /i-quasar SS 433. Thus, at least in this case, the ejecta are made of 
ordinary matter, and not of a fancier substance, such as e + e~ pairs. 

5 The Cannonball Model: summary 

The 'cannon' of the CB model is analogous to the ones of quasars and mi- 
croquasars. In an ordinary core-collapse SN event, due to the parent star's 
rotation, an accretion disk is produced around the newly-born compact object, 
either by stellar material originally close to the surface of the imploding core, or 

by more distant stellar matter falling back after the shock's passage EH) A 
CB made of ordinary-matter plasma is emitted, as in microquasars, when part 
of the accretion disk falls abruptly onto the compact object. Long-duration 
GRBs and non-solar CRs are produced by these jetted CBs. To agree with 
observations, CBs must be launched with LFs, 70^ 10 3 , and baryon numbers 
iV B = O(10 50 ), corresponding to ~ 1/2 of the mass of Mercury, a miserable 
~ 10~ 7 th of a solar mass. Two jets, each with n CB = (n CB ) ~5 CBs, carry 

E icts = 2 n CB 7o N B m p c 2 ~1.5x 10 51 erg, (5) 

comparable to the energy of the SN's non-relativistic shell, that is 0(1%) of 
the explosion's energy, ~98% of which is carried away by thermal neutrinos. 

We have seen that long GRBs are indeed made by SNe, as advocated 
in the CB model well before the pair GRB030329/SN2003dh convinced the 

majority. But do SNe emit cannonballs? Until 2003^3), there was only one 
case with data good enough to tell: SN1987A, the core-collapse SN in the LMC, 
whose neutrino emission was seen. Speckle interferometry data taken 30 and 
38 days after the explosion did show two back-to-back relativistic CBs, see 
fig- [Hl3,f. The approaching one was superluminal: seemingly moving at v > c. 

A summary of the CB model is given in Fig. [T2J The 'inverse ' Compton 
scattering (ICS) of light by electrons within a CB produces a forward beam of 
higher-energy photons: a pulse of a GRB or an XRF. The target light is in a 
temporary reservoir: the glory, illuminated by the early SN light and illustrated 
by analogy in fig. lllH . A second mechanism, synchrotron radiation (SR), takes 
over later and generally dominates the AG. The 7 rays ionize the ISM on the 
CB's path. The CBs collide with the ISM electrons and nuclei, boosting them 
to cosmic ray status. The ISM penetrating a CB's plasma creates turbulent 
magnetic fields within it. The ISM electrons moving in this field emit the 
mentioned SR. This paradigm accounts for all properties of GRBs and CRs. 

The observed properties of a CB's radiation depend crucially on the angle 




Figure 12: The CB model of GRBs, XRFs and CRs. A core-collapse SN 
results in a (black) compact object, a fast-rotating torus of non-ejected material 
and a (yellow) shell of non-relativistic ejecta. Matter (not shown) episodically 
accreting into the central object produces two narrowly collimated beams of 
CBs; some of the 'Northern' CBs are depicted. As the CBs move through the 
glory of non-radial light surrounding the star, they forward Compton up-scatter 
its photons to GRB or XRF energies, depending on how close the line of sight 
is to the CBs' direction. Each CB produces a GRB 'pulse'. Later, a CB gathers 
and scatters ISM particles, which are isotropized by its inner magnetic field. 
In the SN rest system the particles are boosted by the CB's motion: they have 
become CRs. The synchrotron radiation of the gathered electrons is the late 
AG of the GRB or XRF. As the CBs' collisions with the ISM slow them down, 
the CBs generate CRs all along their trajectories, in the galaxy and its halo. 
CRs are also forward-produced, diffusing thereafter in the local magnetic fields. 



6 of its motion relative to the line of sight to the observer, via the Doppler factor 

S = 1/[ 7 (1 - cos 6)} « 27/(1 + V) (6) 

by which a photon's energy is boosted from the CB's rest system to that of a 
(cosmologically nearby) observer. For an isotropic emission in the CB's sys- 
tem, the observed photon number, energy flux and luminosity are ex S, S 3 , S , 
respectively, just as in a v beam from 7r decay. That makes GRBs observable 
only extremely close to one of their bipolar CB axes, 9 = 0(1/7) ~ 1 mrad 
[typically y(t = 0)=7o~ j <5o~10 3 ; and AGs are observed till 7(i)~ 7 o/2]. 

The relation between CB travel-time in the host galaxy, dt+ = dx/{j3c), 
and observer's time, t, is dt+/dt = 'y S/ (1+z). Stop in awe at this gigantic factor: 
a CB whose AG is observed for 1 day may have travelled for 0(1O 6 ) light days, 
what a fast-motion video! A CB with 6*= 1/7= 10~ 3 moves in the sky at an 
apparent transverse velocity of 2000 c, yet another large Doppler aberration. 

6 GRB afterglows in the CB model 

Historically, two GRB phases were distinguished: a prompt one, and the after- 
glow. Swift data have filled the gap, there is no longer a very clear distinction. 
Nor is there a profound difference between the CB-model's radiation mecha- 
nisms, since synchrotron radiation is but Compton scattering on virtual pho- 
tons and, in a universe whose age is finite, all observed photons were virtual. 

In the understanding of GRBs in the CB model, SR-dominated AGs came 
first. The CB-model AG analysis is strictly a 'model': it contains many sim- 
plifications. But the comparison with data determines the distributions of the 
relevant parameters. Given these, the predictions for CRs and for the ICS- 
dominated phase of GRBs (such as all properties of the 7-ray pulses) involve 
only independent observations, basic physics and no 'modeling'. For the reader 
who might want to move to the more decisive sections, I anticipate the con- 
tents of this one. The distribution of 70 and 70 So values of pre-Swift GRBs are 
shown in fig. [T5a,b. The radius of a CB evolves as in fig. [T5H .. A CB does not 
expand inertially; for most of its trajectory it has a slowly changing radius, as 
a common cannonball does. The baryon number of observed CBs is of O(10 50 ). 

To determine the fate of a CB, we make the following assumptions. CBs 
initially expand at (3 S c = C(c/v / 3), the relativistic speed of sound, swiftly be- 
coming spherical in their rest system and losing memory of their initial size. 
For the CB's baryon number returned by the analysis, this means that CBs be- 
come 'collisionless' fast: their nuclei and electrons do not often collide with the 
ISM ones they encounter. Hadron and Thompson cross sections being similar, 
CBs also become transparent, except to long radio waves, losing their radia- 
tion pressure. In agreement with first-principle calculations of the relativistic 
merger of two plasmas ; a chaotic magnetic field is generated within a CB 
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Figure 13: Right to left, top to bottom, a) The distribution of Log 10 (7) for a 
score of pre-Swift GRBs of known z. b) The distribution of a = Log 10 (7<5/10 6 ) 
for the same GRBs, also extracted from the analysis of their AGs. c) The distri- 
bution of peak or 'break' energies in the spectrum of many pre-Swift GRBs. The 
line is the CB-model prediction, based on Eq. (fT2|) and the observed a distribu- 
tion, d) The behaviour of R(j/jo) in two extremes (diffusive or instantaneous) 
for the way the intercepted ISM ions reexit a CB. The coasting behaviour at 
R^R C and final blow-up at 7— *T are of the form i?«i? c (7o//?7) 2 ^ 3 - 




by the ISM particles it sweeps in. In accordance with observations of similar 
plasmas (such as the ISM itself and the CRs it contains) the CB's magnetic 
field is in energy equipartition with the impinging ISM, resulting in: 

B CB (7) = 3 G ( 7 /10 3 ) [n/(lCT 2 cm" 3 )] 1/2 , (7) 

where n is the ISM baryonic number density, normalized to a typical value in 
the 'superbubbles' in which most SNe and GRBs are born. 

In a CB's rest system the motion of its constituents is an inertial memory 
of the initial radial expansion, whose kinetic energy is larger than the one of 
the CB's magnetic field. An ISM proton entering a CB will meander it its 
magnetic field and be isotropically reemitted (in the CB's rest system). The 
rate of radial momentum loss per unit surface is a surface pressure countering 
the expansion. We assume that the dominant effect of this pressure is to 
counteract the expansion. We use Newton's law to compute the ensuing radial 
deceleration and the CB's radius R(j)- The results are shown in fig.fTBTl. A CB 
initially expands quasi-inertially. It subsequently settles into a slowly evolving 
radius till it blows up as its motion becomes non-relativistic (DD06), obeying: 

#(7) w i? c (7o//37) 2/3 , with R c = O(10 14 cm), for typical parameters. (8) 

This is a complex problem, and ours is a big simplification, once assessed by a 
cunning referee as "almost baron Munchhausen" . Yet, the result describes the 
surprising 'jet self- focusing' observed, e.g., in Pictor A, see fig. ITTk. 

The collisions with the ISM continuously decelerate a CB. For a given R(t) 
and ISM baryon number per unit volume n, energy-momentum conservation 
dictates the explicit form of the CB's diminishing Lorentz factor 7(i). Typically 
7(t) is roughly constant for a day or so of observer's time, steepening to oc £ -1 / 4 
thereafter. During the short 7-ray emission time, 7(f) 70. 

We assume that practically all of the energy of the ISM electrons (of 
number density n e ~ n) entering a CB is reemitted fast in the form of SR, 
so that the corresponding observed frequency-integrated AG power per unit 
area is dF / (dt dfi) = n R 2 n e m e c 3 7 2 S 4 / (4 it D 2 l ) , with Dl the luminosity dis- 
tance. The CB deceleration law, dictated by energy-momentum conservation, 
is equally simple: Mq 70 dj = — ir R 2 nm p 7 3 dx, for an H-dominated ISM and 
in the extreme in which the re-emission time of protons is long on the scale of 
the CB's deceleration time. For constant n, the distance travelled by a CB is: 

*(7) « L [( 7 o/7) 2/3 - l] , L = 37V B /(2^i? 2 n 7o 3 ) = 180pc, (9) 

where the number is for N B = 10 50 , n = 10~ 2 cm -3 , R c = 10 14 cm, 70 = 10 3 . 
Given all this, it appears easy to extract from an AG's normalization and 
shape the values of 8, 70 and N s , if one trusts the estimate of i?(7) and uses a 
typical n. Limited observational information makes life a bit harder. 



The spectrum of fig. [7f is actually the one predicted by the CB model, 
illustrated for a typical choice of parameters. The chromatic 'bends' shown as 
dots in this figure, for instance, are 'injection bends': the typical SR energy, in 
the CB's magnetic field, of the electrons entering it, at time t, with a (relative) 
LF 7 = j(t). The small portion of the spectrum above the bend is emitted by 
a tiny fraction of electrons 'Fermi-accelerated' in the CB's turbulent magnetic 
fields to a pre-synchrotron-cooling spectrum E~ p , with 2.2. The predic- 
tion fits with no exception the AGs of the first score of well measured GRBs 
(DDD02/03) of known z. But only on rare occasions can one clearly see in 
an AG the contributions of the various CBs seen in the 7-ray count rate (a 
counterexample is the GRB of figs. [S}d, [5^,b). Thus, generally, the parameters 
extracted from AG fits refer to a dominant CB or to an average over CBs. 

After an observer's day or so, the optical and X-ray AG are typically 
SR-dominated, are above the injection bend, and are of the approximate form: 

F^n^n^Rl^S 4 - 1 ^ 1 - 1 , (10) 

where the unwritten proportionality factors, such as D^ 2 , are known. From 
a fit to the shape of F v (t) one obtains #70 and the combination L of Eq. ©. 
At late times F v oc j(t) 6A cx 7p' 4 with a coefficient determined by the other 
fit parameters. The value of 70 is extracted from the 6.4 root of the inverse 
of Eq. (fT0|) , so that, for a result within a factor of 2, one can tolerate large 
errors in the chosen n or in the estimate of R~. Trusting these, one can extract 
N B from L, perhaps with an uncertainty of one order of magnitudtQ- Eq. (|i0[) 
has been used to fit, after the early fast fall-off, the X-ray and optical data of 
figs. I8I9ITU1 The required form of 7(f) is Eq. ([5]), supplemented by the relation 
between CB's mileage and observer's time, see the end of Sec. 5. 

7 A GRB's 7 rays in the CB model (DD06) 

A pulse of a GRB is made by a CB crossing the parent star's glory. The glory 
is a reservoir of non-radially directed light, fed by the parent star's luminosity, 
as in the artist's view of fig. fTTH . For the best studied GRB-associated SN, 
1998bw, and for 0(ld) after the explosion, the luminosity was L SN ~5xl0 52 
erg/s, in photons of typical energy eV. We adopt these values as 'priors' 

(parameters to be used in calculations, but independent of the CB model). 
Massive stars destined to 'go supernova' eject solar-mass amounts of matter 
in successive explosions during their last few thousand pre-SN years. At the 
'close' distances of O(10 16 cm) relevant here, these stellar coughs generate a 
thick layer of 'wind- fed' material with an approximate density profile pocr~ 2 



1 This is what we did in DDD02/03 but not quite what we wrote. 1 am 
indebted to J. Steinberger for noticing this error. 



and normalization pr 2 ~ 10 16 g cm -1 , the last prior we need. The very early 
UV flash of the SN suffices to ionize the wind-fed matter. The Thompson cross 
section cr T is such that this matter is semitransparent: er T pr 2 /m p — 4x 10 15 
cm. This means that the number of times a SN photon reinteracts on its way 
out -becoming 'non-radial'- is of 0(1), and that the number density of such 
photons is H 7 (r) ~ L SN /(47rr 2 cEi). From emission-time to the time it is still 
one 7-ray interaction length inside the 'wind', a CB has travelled for 

,w_ (Q3c) Pr 2 l + z 10 6 (n) 
*t r -^^l i6 gcm -i 2 l0 5 ( ' 

of observer's time. That is a typical 7-ray pulse rise-time in a GRB, and the 
reason why, closing the loop, distances of O(10 16 cm) were relevant. 

In the collapse of a rotating star, material from 'polar' directions should 
fall more efficiently than from equatorial directions. The CBs would then be 
emitted into relatively empty space. We assume that the wind-material is also 
under-dense in the polar directions. This is not the case for the glory's photons, 
which have been scattered by the wind's matter, and partially isotropised. 
During the production of 7-rays by ICS, 7 — 70- 

Consider an electron, comoving with a CB at 7 = 70, and a photon of 
energy Ei moving at an angle Oi relative to f. They Compton-scatter. The 
outgoing photon is viewed at an angle 9. Its energy is totally determined: 

v5 , , , . l + cos6»j Ei j 6 2 „ , 

E P = X_ d+eos^) Ek = (250 keV) , ^ , s X_ _ , (12) 

where I set w 1 and, for a semi transparent wind, (cosflj) ~ — 1/2. For pre- 
Swift GRBs (z) «1 and, for the typical 7 and S, £' = 250 keV, the average peak 
or 'break' energy in Eq. ©■ From the fits to the AGs of the subset of known 
z, we could determine the distribution of a values, see fig.[T3b. Its fitted result 
is used in Eq. (|T2"]) to predict the overall E p distribution, see fig. [T3T ;. 

The rest of the properties of a GRB's pulse can be derived on similarly 
trivial grounds and with hardly more toil. During the GRB phase a CB is still 
expanding inertially at a speed (3 S c. It becomes transparent when its radius is 
i?tr ~ [3 ct t A^ B /(47r)] 1 / 2 , at an observer's time very close to that of Eq. pT|) . for 
typical parameters. One can simply count the number of ICS interactions of 
a CB's electrons with the glory, multiply by their energy, Eq. (fl"2]) . and figure 
out the isotropic-equivalent energy deduced by an observer at an angle 0: 



g^^VM >^^3. 2xl0 53erg, (is) 
6 c V 4i 



where the number is for our typical parameters, and agrees with observation. 



7.1 ^ Is it Inverse Compton Scattering ... 

The 7 and 6 dependance in Eqs. (|12I13|) is purely 'kinematicaP, but specific to 
ICS: it would be different for self-Compton or synchrotron radiation. To verify 

that the 7 rays of a GRB are made by ICS, as proposed by Shaviv and 

Dar, we may look at the correlations between GRB observables. 

In the CB model, the (7, 5, z) dependences of the peak isotropic luminosity 
of a GRB, Lp°; its pulse rise-time, t r i SC ; and the lag-time between the peaks of 

a pulse at different energies, ti ag ; are also simply derived to be: 

(1 + z) 2 L p so oc S\ t lise oc(l + z)/( 7 S), t lag cc(l + z) 2 /(<5 2 7 2 )- (14) 

I have not specified the numerical coefficients in Eqs. (|14[) . which are explicit, as 
in Eqs. (|12I13[) . Of all the parameters and priors in these expressions, the one 
explicitly varying by orders of magnitude by simply changing the observer's an- 
gle is (5(7, 9), making it the prime putative cause of case-by-case variability. For 
such a cause, Eqs. {T3J) and the first of QU imply that £^ so oc [(1 + zf Ljf°] 3 / 4 . 
This is tested in fig. [7^. A most celebrated correlation is the [E p , -E^ so ] one, see 
Fig.[7p. It evolves from E p oc [E'^ } 1 / 3 for small E p , to E p oc [£^ so ] 2/3 for large 
E p . This is because the angle subtended by a moving CB from its place of 
origin is (3 s /j, comparable to the beaming aperture, I/7, of the radiation from 
a point on its surface. Integration over this surface implies that, for #<C I/7, 
S oc 7, while in the opposite case S varies independently. The straight lines in 
fig. \7jp> are the central expectations of Eqs. (|12ll3p . the data are fit to the pre- 
dicted evolving power law. The predicted [ti ag , L p so ] and [t T i Se , L p so ] correlations 
are tested in figs. [7J:,d. The seal of authenticity of inverse Compton scattering 
— by a quasi-point-like electron beam — is unmistakable in all of this, QED. 

7.2 ... on a Glory's light ? 

The 'target' photons subject to ICS by the CB's electrons have very specific 
properties. Their number-density, n 7 (r) oc L SN /r 2 , translates into the ~ t~ 2 
late-time dependence of the number of photons in a pulse since, once a CB is 
transparent to radiation, ICS by its electrons simply 'reads' the target-photon 
distribution. As a CB exits the wind-fed domain, the photons it scatters are 
becoming more radial, so that l + cos#i — > r -2 oc t~ 2 in Eq. (fT2")l . For a semi- 
transparent wind material, which we have studied in analytical approximations 
and via Montecarlo, this asymptotic behaviour is reached fast and is approxi- 
mately correct at all t. This means that the energies of the scattered photons 
evolve with observer's time as t~ 2 : the 'Et 2 law' of Eq. ^ and fig. [7^. 



7.3 The pulse shape and the spectrum 

The spectrum of a GRB, Eq. ([3]), and the time-dependence of its pulses, Eq. ((4]), 
describe the data well, and are actually analytical approximations to the results 
of ICS of an average CB on a typical glory. The spectrum of a semitransparent 
glory has a 'thermal bremsstrahlung' shape, dn 7 /dEi(x (TijEi) a Exp[— Ei/Ti], 
with a ~ 1 and Ti ~ 1 eV. The first term in Eq. ([3]) is this same spectrum, 
boosted by ICS as in Eq. (fT2|) . by electrons comoving with the CB, E e =-fm e c 2 . 
The second term is due to ICS by 'knock-on' electrons (generated while the CB 
is not yet collisionless) and electrons 'Fermi-accelerated' by the CB's turbulent 
magnetic fields. They both have a spectrum dn e /dE e (xE~P, with /3~2 to 2.2. 
They are a small fraction of the CB's electrons, reflected in the parameter b, 
which we cannot predict. The temporal shape of a pulse has an exponential rise 
due to the CB and the windy material becoming transparent at a time ~ £J , 
see Eq. (JTTJ), the width of pulse (in 7 rays) is a few tf T , the subsequent decay 
is oct -2 . The time-energy correlations obey the l Et 2 law'. All as observed. 

7.4 Polarization 

A tell-tale signature of ICS is the high degree of polarization. For a pointlikc 
CB the prediction^ is n«2 6» 2 7 2 /(l+6» 4 7 4 ), peaking at 100% at 6 = l/ 7 ,thc 
most probable 9, corresponding to 90° in the CB's system. For an expanding 
CB, n is a little smaller. For SR, which dominates the AGs at sufficiently late 
times, the expectation is II w 0. The 7-ray polarization has been measured, 
with considerable toil, in 4 GRBs. It is always compatible, within very large 
errors, with 100%. The situation is unresolved^!'. I shall not discuss it. 

8 Detailed Swift light curves and hardness ratios 

Swift has abundantly filled its goal to provide X-ray, UV and optical data 
starting briefly after the detection of a GRB: compare the Swift result of fig. 8a 
to the pre-Swift data in fig. 8b. In the CB-model description of the data in 
figs. 8,9,10, the abruptly falling signal is the tail of one or several 7-ray pulses or 
X-ray flares, produced by ICS and jointly described by Eqs. (|3I4[) . The following 
'afterglow', its less pronounced decay and subsequent achromatic 'bend' are due 
to the CBs' synchrotron radiation, described by Eq. (fT0|) . Thanks to the quality 
of SWIFT data one can proceed to test these CB-model predictions in detail. 

The two prompt optical 'humps' of GRB 060206 in fig. [5^ are the ICS 
low-energy counterparts of its two late X-ray flares of fig. [5Jd, simultaneously 
fit by Eqs. (|3I4|) . Swift provides a rough measure of a GRB's spectrum: the 
hardness ratio of count rates in the [1.5-10] keV and [0.3-1.5] keV intervals. 
Given the case-by-case parameters of a CB-model fit to the [0.3-10] keV light 



curve, one can estimate the corresponding hardness ratio ^* . This is done 
in figs. [9h,d and [10ji,c for GRB 060904 and XRF 060218, respectively. This 
last XFR is observed at a 'large' angle, 9 ~ 5 mrad and a correspondingly 
small So, its single X-ray pulse is, in accordance with Eq. (fTTj) . relatively wide. 
The optical and UV counterparts of the X-ray pulse are clearly visible as the 
'humps' in the optical data of fig. [T0b . Given the 'Et 2 law' of Section 7.3, 
the pulse peak times at different frequencies are simply related: t pca kOc_E -1 / 2 . 
The prediction, an example of the ubiquitous 1/r 2 law of 3-D physics, is tested 

in fig. [7^. The peak fluxes at all frequencies are also related as dictated ® by 
Eq. Q . The adequacy of the CB model over many decades in flux and time is 
exemplified by the X-ray light curve of GRB 061121 in fig. [TOH. 

The predictions for the peak 7-ray energy of Eq. (|T2|l , its distribution as 
in fig. [T3"b. the GRB spectrum of Eq. ([3]), and the correlations of figs.[7^-d are 
clear proof that ICS is the prompt GRB mechanism. The test of the Et 2 law in 
fig- 03 corroborates that the 'target light' becomes increasingly radially directed 
with distance: Inverse Compton Scattering on a 'glory's light' by the electrons 
in CBs is responsible for the j-ray pulses of a GRB and their sister X-ray flares 
and optical humps. The properties of the subsequent synchrotron-dominated 
afterglows are also in accordance with the CB model. 

9 The GRB/SN association in the CB model 

We have gathered very considerable evidence that the LFs and viewing angles 
of observed GRBs are 7o =O(10 3 ) and = 0(1) mrad. The fraction of GRBs 
beamed towards us is ~ 9 2 = O(10~ 6 ). The number of such observed GRBs 
(with a hypothetical Air coverage) is a few a day. The same coverage would 
result in the observation of a few million core-collapse SN per day, in the 
visible Universe. These numbers are compatible with the extreme conclusion 
that all these SNe emit GRBs, but the estimates and errors are sufficient 
to accommodate a one order of magnitude smaller fraction, which would be 
compatible with most Type Ib/c emitting (long) GRBs. 

10 Short Hard 7-ray Bursts (SHBs) 

SHBs share with (long) GRBs the properties not reflected in their name. A 
good fraction of SHBs have 'canonical' X-ray light curves. The origin of SHBs 
is not well established, in contrast to that of GRBs and XRFs. Clues to the 
origin and production mechanism of SHBs are provided by their similarity to 
long GRBs. The X-ray light curves of some well-sampled SHBs are 'canonical'. 
The similarities suggest common mechanisms generating the GRB and SHB 
radiations. This is expected in the CB model, wherein both burst types are 
produced by highly relativistic, narrowly collimated, bipolar jets of CBs, ejected 



in stellar processes . The mechanisms for their prompt and AG emissions 
(ICS and synchrotron) coincide with the ones of GRBs. The 'engine' is different; 
it is a core-collapse supernova for GRBs and XRFs, in SHBs it may be a merger 
(of two neutron stars or a neutron star and a black hole), the result of mass 
accretion episodes on compact objects in close binaries (e.g. microquasars), or 
phase transitions in increasingly compactified stars (neutron stars, hyper-stars 
or quark stars), induced by accretion, cooling, or angular- momentum loss. 

In the CB model, the 'master formulae' describing prompt and afterglow 
emissions in long GRBs are directly applicable to SHBs, provided the parame- 
ters of the CBs, of the glory, and of the circumburst environment, are replaced 
by those adequate for SHBs. This results in a good description of the data^P. 

11 Cosmic Rays in the CB model 

In the CB model, CRs are as simple to understand as GRBs. If relativistic 
CBs are indeed ejected by a good fraction of core-collapse SNe, it is inevitable 
to ask what they do as they travel in the ISM. The answer is that they make 
CRs with the observed properties, simply by interacting with the constituents 
of the ISM, previously ionized by the 7 rays of the accompanying GRB. Early 
in their voyage, CBs act as Compton relativistic rackets, in boosting a glory's 
photon to 7-ray status. Analogously, all along their trajectories, CBs act as 
Lorentz relativistic rackets, in boosting an ISM nucleus or electron to CR status. 
Once again, the necessary input is two-fold. On the one hand, there are the 
properties of CBs: the average number of significant GRB pulses (or CBs) per 
jet (5), the 70 distribution of fig. [T3"a , and the N B ~ 10 50 estimate. On the 
other hand, there are a few 'priors', items of information independent of the 
CR properties: the rate of core-collapse SNe, the relative abundances, n A (of 
the elements of atomic number A) in the ISM, and the properties of Galactic 
magnetic fields. 

We shall see that the CB-model predictions for the normalization of CR 
spectra are correct to within a factor of 0(3), while the ratios between elements 
are correct within errors. In figs. 1, 3 and 4a, the predictions have been made to 
adjust the data, not reflecting the common overall normalization uncertainty. 

11.1 Relativistic rackets: The knees 

Our simplest result concerns the 'knees' of the all-particle spectrum in fig. |TJl 
and of the main individual elements in fig. [3] The essence of their understanding 
is kinematical and trivial. In an elastic interaction of a CB at rest with ISM 
electrons or ions of LF 7, the light recoiling particles (of mass mw Am p ) retain 
their incoming energy. Viewed in the ISM rest system, they have, for large 
7, a flat spectrum extending up to E ~ 2j 2 mc 2 [this is recognizable as the 



forward, massive-particle, z = 0, analog of Eq. (fT2"|) ]. Thus, a moving CB is 
a gorgeous Lorentz-boost accelerator: the particles it elastically scatters reach 
up to, for 7 = 70 = (1 to 1.5) x 10 3 , an A-dependent knee energy Ek nee (A) ~ 
(2to4)x 10 15 AeV. If this trivial process is the main accelerator of CRs, there 
must be a feature in the CR spectra: endpoints at E^e^A). The arrows in 
fig. [3] show that the H and He data are compatible with this prediction. So 
does the second knee of fig. QJi, the predicted Fe knee. The CR flux above the 
H knee, to which we shall return, is ^10 -15 of the total. 

11.2 The spectra below the knee 

The 'elastic' scattering we just described is dominant below the knees. To com- 
pute the resulting spectrum, we assume that the ISM particles a CB intercepts, 
trapped in its magnetic mesh, reexit it by diffusion, isotropically in the CB's 
system, and with the same 'confinement' law, Eq. |T]), as in the Galaxy (the 
opposite assumption, that they are immediately elastically scattered, yields a 
slightly different spectral index). The CB deceleration law is Eq. its radius 
evolves as in Eq. ((SJ. A modest amount of algebra gives a simple result (DD06), 
which, for 7 > 2 and to a good approximation, rea 

d 7ca A \ Z ) Jl 7 7/3 7max[7,7 CR /(2 7)] 7co 

where 7 CR is the CR's LF, and (3 con { is the same confinement index as in 
Eq. ([1]). The flux dF e i ast /dj CR depends on the priors n A , /3 ccm i, and 70, but 
not on any parameter specific to the mechanism of CR acceleration. But for 
the normalization, this flux is A-independent. In the large range in which it is 
roughly a power law, cLF elast /d7 CR oc [7 CR ] _/3cR , with /3 CR =13/6«2.17. 

The H, He and Fe fluxes of fig. [3] are given by Eq. (fT5")) . modified by the 
Galactic confinement r-dependence of Eq. (fT]), with /3 con f = 0.6. The fastest- 
dropping curve in fig.[3£i corresponds to a fixed 70. The other two curves are for 
the 70 distribution of fig. 113b . and one twice as wide. The low-energy data of 
fig-lb are also described by Eq. (fT5"]) . whose shape in this region (the 'hip', also 
visible in fig.[3j: for Fe) is insensitive to 70 and, thus, parameter-independent. 

11.3 The relative abundances 

It is customary to discuss the composition of CRs at a fixed energy E A = 1 
TeV. This energy is relativistic, below the corresponding knees for all A, and 
in the domain wherein the fluxes are dominantly elastic and well approximated 

2 I am keeping factors of A/Z for kicks. Numerically, they are irrelevant: 
the theory and data are not so precise, and (A/Z) - 6 is 1 for H, 1.6 for Fe. 



by a power law of index (3th = /3oiast+/3conf — 2.77. Expressed in terms of energy 
(E A (xAj), and modified by confinement as in Eq. (Q]), Eq. (fT5|) becomes: 

dF ohs /dE A <x n A A^- 1 E~^\ X CR (A) = (njn p ) A 1 - 77 , (16) 

with n A an average ISM abundance and X CR (A) the CR abundances relative 
to H, at fixed E. The results, for input n A s in the 'superbubbles' wherein most 
SNe occur, are shown in Fig. [2Jd. In these regions, the abundances are a factor 
~ 3 more 'metallic' than solar (a 'metal' is anything with Z > 2). Eq. (| 16[) snugly 
reproduces the large enhancements of the heavy-CR relative abundances, in 
comparison with solar or superbubble abundances (e.g. A 1 - 77 = 1242 for Fe). 
The essence of this result is deceptively simple: in the kinematics of the collision 
of a heavy object (a CB) and a light one (the ISM nucleus), their mass ratio 
(N B /A ^oo!) is irrelevant. 

11.4 Above the knees 

We discussed around Eq. ([7]) the generation of turbulently-moving magnetic 
fields (MFs) in the merger of two plasmas. Charged particles interacting with 
these fields tend to gain energy: a relativistic-injection, 'Fermi' acceleration 
process, for which numerical analyses^® result in a spectrum dN/dEocE~ 2 - 2 , 
p~ 2.2. For the ISM/CB merger, we (DD04) approximate the spectrum of 
particles accelerated within a CB, in its rest system, as: 

dN/d lA ex 7 ; 2 ' 2 6( 7a - 7 ) e[ 7roax ( 7 )- 7 J, 7max ~ 10 5 7o 2/3 (Z/A) y/ 3 , (17) 

The first function reflects the fact that it is much more likely for the light 
particles to gain than to lose energy in their elastic collisions with the heavy 
'particles' (the CB's turbulent MF domains). The second is the Larmor 
cutoff implied by the finite radius and MF of a CB, with a numerical value 
given for the typical adopted parameters. But for the small dependence of 7max 
on the nuclear identity (the factor Z/A), the spectrum of Eq. (fl?) is universal. 
Boosted by the CB's motion, an accelerated and re-emitted particle may reach 
a Larmor-limited 7ca [max] =2 77max , a bit larger, for 7 = 7 o~ 1.5 x 10 3 , than 
the corresponding GZK cutoffs. Our model has a single source, CBs, for the 
acceleration of CRs from the lowest to the largest observed energies. 

The calculation of the 'elastic' spectrum of Eq. (fT5)) was done for the bulk 
of the ISM particles entering a CB, assuming that they were not significantly 
Fermi-accelerated, but kept their incoming energy, i.e. dN/dj A oc <5( 7a — 7 ). 
The 'inelastic' spectrum, with dN/d r y A as in Eq. (|15p . yields an equally simple 
result. The two i? 2 -weighed spectra are shown (for H) in fig. HU The inelastic 
contribution is a tiny fraction of the total, and is negligible below the knee, a 
point at which we may compare the ratio of fluxes, /, the only parameter freely 
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Figure 14: Contributions to the i5 2 -weighed proton source spectrum. 

fit to the CR data. The boost of ISM particles by a CB and their acceleration 
within it are mass-independent, so that the ratio / is universal. 

The i? 3 -weighed source spectra for the main elements are shown in fig. 115b . 
They are very different from the observed spectra of fig. 115b . for many reasons. 
Below the ankle(s) the slopes differ due to Galactic confinement, see Eq. ((T|). 
Above the ankles the flux from Galactic sources is strongly suppressed: we 
would see their straight-moving CRs only for CB jets pointing to us. The 
CRs above the ankle are mainly extragalactic in origin, and they also cross the 
Galaxy just once. Extragalactic CRs of A > 1 arc efficiently photo-dissociated 
by the cosmic infrared light. Extragalactic CRs are GZK-cutoff. All this can 
be modeled with patience and fair confidence. Below the ankle extragalactic 
CRs may have to fight the CR 'wind' of our Galaxy, analogous to that of the 
Sun. We have covered our lack of information on this subject by choosing two 
extreme possibilities (DD06), resulting in the two curves of figs.[lji and[3Ji,b. 

In fig. QJ) I have converted the results of fig. [15b into a prediction for 
(LiiA(E)). The flux at the second knee is dominated by Galactic Fe at its 
knee. Thereafter this flux decreases abruptly to let extragalactic H dominate 
all the way from the ankle to the nominal position of the proton's GZK cutoff. 
Above that point the high-energy tail of Galactic Fe may dominate again. 



11.5 The CR luminosity, and the overall normalization of the CR flux 

The rate of core-collapse SNe in our Galaxy is i? SN ~ 2 per century. In the CB 
model, we contend that ~ 50% of the energy of CRs is transfered to the mag- 
netic fields they generate . If all core-collapse SNe emit CBs, the Galactic 
CR luminosity should be L CR ~ i? SN Ej cts /2 ~ 4.7 x 10 41 erg s _1 , with Ej cts as 
in Eq. ((SJ). This is 3 times larger than the rhs of Eq. @. The 'discrepancy' is 
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Figure 15: Predicted spectra for the abundant elements and groups. The 
vertical scales are E 3 dF/ dE. (a): The source spectra, with a common arbitrary 
normalization, (b): The normalized CR spectra at the location of the Earth. 
Notice that the horizontal and vertical scales are different in (a) and (b). 



not worrisome. A smaller fraction of SNe may generate high-70 CBs. The rhs 
of Eq. © is for 'standard' CRs, but the confinement volume and time of the 
CB model are non-standard by factors of ~10. All inputs are fairly uncertain. 

The calculation of the flux above the ankle is lengthy but straightforward. 
But for the GZK effect, its shape is that of the source H flux, since protons at 
that energy should escape other galaxies directly, and enter ours unhindered. 
Its normalization per SN is fixed. The SN rate per unit volume is measured in 
the local Universe. The overall flux is the result of the integration over redshift 
of the flux from past SNe. The integrand must be properly red-shifted and 
weighed with the star-formation rate as a function of z (SN progenitors have 
short lives on Hubble-time scales). The integration in z is an integration over 
look-back time, as opposed to distance, since CRs do not travel straight. The 
error in the result is hard to estimate, its central value is within a factor of 2 of 
the observations (DD06). This explains the coincidence that the ankle is the 
escape energy of protons from the Galaxy and the place where the extragalactic 
flux -not enhanced by confinement and thus less steep- begins to dominate. 

11.6 CR diffusion, CR electrons, and the 7 background radiation 

In the standard paradigm CRs are accelerated by the nonrelativistic ejecta of 
SNe. SNe occur mainly in the central realms of the Galaxy, so that CRs must 
diffuse to arrive to our location. A directional asymmetry is predicted, and not 

observed . For CR electrons the problem is even more severe: their cooling 
time in the Galaxy's light, and magnetic fields, is so short that they should 
have lost all their energy on their way here. 

Our source distribution is totally non-standard, CBs generate CRs all 
along their many-kpc-long trajectories, see Eq. © and fig. [T2"b . Depending on 
the ISM density profile they encounter, CBs may travel for up to tens of kpc 
before they become nonrelativistic. It takes some 6xl0 4 years to travel 20 kpc 
at u~c. If a Galactic SN occurs every 50 years, and emits an average of 10 
CBs, there are currently several thousand CBs in the Galaxy and its halo. This 
is a very diffuse CR source, satisfactory in view of the previous paragraph. We 
have not yet studied the CR source-distribution and diffusion in detail. 

Below their knee at 2jQjn e c 2 ~ 2.3 TeV, the source spectrum of CB- 
accelerated electrons has the same index as that of nuclei: dN/dE e oc E~"cr , 

/3 ca « 2.17. The predicted spectrum ^3, steepened by radiative energy loses, 
has an index (3 e = /3 CR +1~ 3.17. Its observed slope ^3) [ s 3.2 ± 0.1 above 
.E e ~10 GeV, an energy below which other losses should dominate (DD06). 

The Gamma Background Radiation (GBR), measured by EGRET from 
a few MeV to ~ 10 5 MeV, was argued to be dominantly of cosmological origin, 

in directions above the disk of the Galaxy and away from its bulge . A 
more careful analysis reveals a significant correlation of its intensity with our 



position relative to the Galactic centre ^3 . The CB-model reproduces this 
correlation, provided a good fraction of the GBR is generated by CR electrons 
at high galactic latitudes, as they cool radiatively by the very same process that 
steepens their spectrum. The predicted index of the radiated GBR photons is 

/ a y = (/? e -l)/2 = 2.08. The observed one^U) is 2.10±0.03. 

12 If CBs are so pervasive, why are they not readily observed? 

The answer is simple. The CBs of SNe are tiny astrophysical objects: their 
typical mass is half of the mass of Mercury. Their energy flux at all frequencies 
is cx<5 3 , large only when their Lorentz factors are large. But then, the radiation 
is also extraordinarily collimated, it can only be seen nearly on-axis. Typically, 
observed SNe are too far to photograph their CBs with sufficient resolution. 

Only in two SN explosions that took place close enough, the CBs were 
in practice observable. One case was SN1987A, located in the LMC, whose 
approaching and receding CBs were photographed, see fig. [5^,f. The other case 
was SN2003dh, associated with GRB030329, at z = 0.1685. In the CB model 
interpretation, its two approaching CBs were first 'seen', and fit, as the two- 
peak 7-ray light curve of fig. [SJd and the two-shoulder AG of fig. [5^,b. This 
allowed us to estimate the time-varying angle of their superluminal motion in 
the sky 113). Two sources or 'components' were indeed clearly seen in radio 
observations at a certain date, coincident with an optical AG rebrightcning. 
We claim that the data agree with our expectation^, including the predicted 

inter-CB separation 03 of fig. [BJi. The observers claimed the contrary, though 
the evidence for the weaker 'second component' is > 20a. They report 
that this component is 'not expected in the standard model'. The unpublished 
and no-doubt spectacular-discovery picture of the two superluminally moving 
sources would have been worth a thousand words... in support of the CB model. 

13 Other CR sources 

We have defended the simplistic view that CBs from SNe is all one needs to 
generate CRs at all energies. The recent data of Auger ED show a signifi- 
cant correlation between the arrival directions of ultra-high energy CRs (UHE- 
CRs) and Active Galactic Nuclei (AGNs) located within a distance of 75 Mpc 

3 The size of a CB is small enough to expect its radio image to scintillate, 
arguably more than observed EH) , Admittedly, we only realized a posteriori 
that the ISM electrons a CB scatters, synchrotron-radiating in the ambient 
magnetic field, would significantly contribute at radio frequencies, somewhat 
blurring the CBs' radio image ^D. 



(2 ~ 0.02), comparable to the GZK 'horizon' at the observed energies, E > 56 
EeV. As the authors discuss, this does not mean that AGNs are the actual 

sources , for AGNs are themselves correlated with matter and with active 
regions of enhanced stellar birth and death. More work on correlations is no 
doubt in progress. A search for correlations with GRBs is less hopeful, for the 
fraction of them observed from within z ~ 0.02 is negligible, and the observed 
'correlated' CRs may have been bent by magnetic fields up to a few degrees, 
implying a possible delay of millenia between the arrival times of 7's and CRs. 

Observations of AGNs are an ingredient of the 'inspiration' of the CB 
model, as we have discussed in connection with Pictor A, see fig. [TTa.b. Nat- 
urally, we have estimated their contribution to the UHECR proton flux, con- 
cluding that they may constitute at most 1 to 10% of the flux generated by 
extragalactic SNe (DD04). The estimate is for the energy-integrated flux; in 
applying it to UHECRs we assumed the same energy dependence for the flux 
generated -by the same mechanisms- by the CBs of SNe and AGNs (a small 
difference of spectral index implies an enormous uncertainty). In a subsequent 

study ^S) f the CR electron flux, assumed to be in a fixed proportion to the 
proton flux, we used more recent inputs, and simplified and modified our upper 
limit to 40%. But we forgolQ to extract the putative consequences for the AGN 
contribution to UHECR protons! 

14 Discussion and conclusion 

We do not have a solid understanding of accretion onto black holes or neutron 
stars. But such processes are observed to result in the ejection of relativistic 
and highly collimated jets. We assumed that a similar process takes place as a 
stellar core collapses, leading to a supernova event. We posited that the SN's 
relativistic ejecta -two jets of cannonballs- are the sources of GRBs and CRs. 
The association between SNe and (long) GRBs is now established. We argued 
that the electrons in a CB, by inverse Compton scattering on the illuminated 
surroundings of the exploding star, generate the 7 rays, X-rays, UV and optical 
light of the 'prompt' phase of a GRB. The ensuing results are the simplest and 
most predictive, they are a firm 'theory'. In this paper I have, however, followed 
the historical development, in which the CB parameters were first extracted 
from the observations of the afterglow of GRBs. This involves a 'model', a set 
of arguable but simple hypothesis leading to the prediction of the properties 
of the AG -dominated by synchrotron radiation by the ISM electrons that a 
CB intercepts- as a function of frequency and time. In the historical order 

4 In these days of large experimental collaborations percolated by theorists, 
rumours herald publications. It might have been possible to turn this comment 
into a renewed timely 'prediction', prior to the Auger announcement. 



the 'prompt' results for GRBs are predictions of the theory. Some results for 
Cosmic Rays -the ISM particles that CBs scatter in their journey- are also 
'theory', others can be viewed as further tests of the 'model'. 

The results for GRB afterglows may be based on a simplified model, but 
they work with no exception all the way from radio to X-rays (DDD02,03). In 
particular, they describe correctly GRB 980425, located at a redshift two orders 
of magnitude closer than average. Its associated SN is the one we 'transported' 
to conclude -thanks to the reliability of our AG model- that core-collapse SNe 
generate long GRBs (DDD02). The X-ray light curve of GRB 980425 and a few 
others, with extremely scarce data, was described with the 'canonical' proper- 
ties later observed in detail in many SWIFT-era GRBs. It is not recognized 
that the two CBs of GRB 030329 were seen, or that their separation in the sky 
was the predicted 'hyper-luminaP one. In view of the overall success of the CB 
model, this is a durable hurdle: GRBs so close and luminous are very rare. 

The accuracy of the predictions for the prompt phase of GRBs amazes 
even the CB-model's proponents. The typical values and the correlations be- 
tween the 7-ray prompt observables leave little doubt that the production mech- 
anism is inverse Compton scattering on 'ambient' light of ~ 1 eV energy. The 
approximate scaling law E dN 1 / dEdt otF(Et 2 ) -spectacularly confirmed in the 
case of XRF 060218- demonstrates that the light is that of a 'glory': the early 
SN light scattered by the 'windy' pre-SN ejecta. A GRB spectrum that works 
even better than the phenomenological 'Band' expression is also predicted. 
The flux and its spectral evolution during the prompt and rapid-decline phases 
are the expected ones, as we tested in minute detail with SWIFT data. 

In the internal-external fireball model of GRBs, highly relativistic thin 
conical shells of e + e~ pairs, sprinkled with a finely tuned baryon 'load', col- 
lide with each other generating a shock that accelerates their constituents and 
creates magnetic fields. Each collision of two shells produces a GRB pulse by 
synchrotron radiation. The ensemble of shells collides with the ISM to produce 
the AG by the same mechanism. The energy available to produce the GRB 
pulse -as two shells moving in the same direction collide- is more than one 
order of magnitude smaller than that of the merged shells as they collide with 
the ISM at rest. The ratio of observed GRB and AG energies is more than one 
order of magnitude, but in the opposite direction. This 'energy crisis' in the 
comparison of bolo metric prompt and AG fluences © is not resolved. More- 
over, the GRB spectrum cannot be accommodated on grounds of synchrotron 
radiation ^HJ) ; the- 'standard' prompt mechanism. The SWIFT-era observations 
also pose decisive problems to the standard model, whose microphysicsl^Sl) ; re- 
liance on shocks anc i correlations based on the jet-opening angle ^^HJ) have 
to be abandoned, according to the cited authors. 

In spite of the above, the defenders of the fireball model are not discour- 



aged. Their attitude towards the CB model, whose observational support is so 
remarkable, is not equally supportive. This may be due to cultural differences. 
Particle physicists believe that complex phenomena may have particularly sim- 
ple explanations. They thrive on challenging their standard views. Doubting 
or abandoning a previous consensus in astrophysics is less easy. 

The CB-model description of Cosmic Rays is also simplistic: there is only 
one source of (non-solar) CRs at all energies, and only one parameter to be 
fit. The model has a certain inevitability: if CBs with the properties deduced 
from GRB physics are a reality, what do they do as they scatter the particles 
of the interstellar medium? We have argued that they transmogrify them into 
CRs with all of their observed properties. The mechanism is entirely analogous 
to the ICS responsible for a GRB's prompt radiation. Suffice it to substitute 
the CB's electron, plus the ambient photon, by a moving-CB's inner magnetic 
field, plus an ambient nucleus or electron. 

After a century of CR measurements, the CB-model results lack the glam- 
our of predictions. Yet, the expectations for the knee energies, and for the 
relative abundances of CRs, are 'kinematical', simple, and verified. They con- 
stitute evidence, in my opinion, that the underlying model is basically correct. 
The prediction for the shape of the spectra: the low energy hips, the large 
energy stretch very well described by a power-law of (source) index /3 S = 13/6, 
and the steepening at the knees, are also verified. The index j3 a is measured 
well enough for the adequacy of the prediction to be sensitive to the details of 
the underlying model, such as the form of the function R("f) in Eq. (J5J) . I can- 
not claim that the fact that the prediction is right on the mark is much more 
than a consistency test, for the physics underlying this aspect of the problem 
may be terrifyingly complex. The same CR source -cannonballs from super- 
novae, this time extragalactic- satisfactorily describes the CR data above the 
ankle. Finally, the properties of CR electrons, and of the high-latitude Gamma 
'Background' Radiation, are also correctly reproduced. 

Most CR scholars agree with the 'standard' paradigm that the flux well 
below the knee is produced by the acceleration of the ISM in the frontal shocks 
of the nonrelativistic ejecta of SNe. In spite of recent observations of large 
magnetic fields ^3) i n collisions of SN shells and molecular clouds, nobody has 
been able to argue convincingly that this process can accelerate particles up to 
the (modest) energy of the knee, and to show that the number and efficiency 
of the putative sources suffices to generate the observed CR luminosity (to my 
satisfaction, I add, to make these statements indisputable). From this point on, 
there is no 'standard' consensus on the origin of CRs, e.g., of the highest-energy 
ones. In this sense, the CB model is regarded as yet another model, which it is. 
After all, we are only saying that CRs are accelerated by the jetted relativistic 
ejecta of SNe, as opposed to the quasi-spherical, non-relativistic ones. Yet, 

the CB model is also rejected by the CR experts, sometimes even in print©, 



though it survives the critique SOJ) But, concerning CRs, the model does not 
trigger the same indignant wrath as in the GRB realm. 

I have shown that the problem of GRBs is convincingly -i.e. predictively- 
solved and that, on the same simple basis, all properties of CRs can be easily 
derived. Only an overwhelmed minority recognizes these facts, in contradic- 
tion with Popper's and Ockham's teachings. I would conclude with a dictum 
attributed to Lev Landau: 'In astrophysics, theories never die, only people do. ' 
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